Topological Chern indices in molecular spectra 
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Topological Chern indices are related to the number of 
rotational states in each molecular vibrational band. Modifi- 
cation of the indices is associated to the appearance of "band 
degeneracies" , and exchange of rotational states between two 
consecutive bands. The topological dynamical origin of these 
indices is proven through a semi-classical approach, and their 
values are computed in two examples. The relation with the 
integer quantum Hall effect is briefly discussed. 

PACS numbers: 03.65. Sq; 02.40.-k; 31.15.Gy; 



Topological numbers play an important role in many 
area of physics B, but their appearance in molecular 
physics and especially in rovibrational problems has not 
been systematically appreciated so far. Simple molec- 
ular systems typically allow adiabatic separation of vi- 
brational and rotational motion. For non-degenerate iso- 
lated electronic state (this is the case of ground state 
for most molecules) the rovibrational energy level system 
consists of vibrational bands, each associated with one or 
several degenerate vibrational states. If the rovibrational 
coupling is not too strong, further splitting of the rovi- 
brational structure into sub-bands can be clearly seen. 
The well-known example is the splitting of the triply de- 
generate vibrational structure for a spherically symmet- 
rical molecule into three sub-bands due to the first-order 
Coriolis interaction Within each sub-band formed 

by 2j + 1 — C levels, with respectively C — +2, 0, —2, all 
energy levels are usually characterized by the quantum 
number j of the total angular momentum, and by another 
quantum number R = j + C/2 which characterizes the 
coupling of j with the vibrational angular momentum. 

In this letter we show that the integer C can be de- 
fined in much more general situation as an additional 
quantum number having a precise topological meaning, 
namely a Chern index, whose construction will be ex- 
plained below. This index is defined in the classical limit 
of the rotational motion. It can be associated with any 
vibrational band presented in the energy level pattern of 
molecular systems. Theorem (||) relates this topological 
index to the number of rotational states within the band. 
A modification of the index is associated with the forma- 
tion of a contact (a degeneracy) between two consecutive 
vibrational bands, and is shown to generically imply an 
exchange of one rotational state between the two bands. 



Such a relation was first conjectured in 1988 Q af- 
ter the study of the simple model (jl]) , and a number of 
effective Hamiltonians reconstructed from experimental 
data (see Ref. ||-f7j for the molecular examples: SiH4, 
CD 4 , SnH 4 , CF 4 , Mo(CO) 6 ). The universal character of 
the redistribution phenomena and its relevance to integer 
Hall effect was discussed on several occasions |§Hl0| . 

To explain the physical phenomenon and to prepare 
the formulation of a rigorous statement let us consider 
a toy problem which involves two quantum angular mo- 
menta J and S, with fixed modulus J 2 = j(j + 1) and 
S 2 = s(s+l) with j, s integer or half-integer: J acts in the 
space Tij of dimension (2j + 1) which is the irreducible 
representation space of the SU (2) group. Similarly, S 
acts in the space 7i a of dimension (2s + 1). The total 
space is Wtc-t — 'Hj^'Hs with dimension (2j + l)(2s + 1). 

The most general quantum Hamiltonian H(S,3/j) we 
will consider is a hermitian operator acting in 'H-iqI and 
its action in space Hj is supposed to be expressed in 
terms of the operators 3/j. The factor 1/j is introduced 
here to ensure the existence of the classical limit for j — > 
oo. An extremely simple form of H is 



H=(l-t)S z +-(J-S), 
J 



(1) 



with j > s, t £ M, which was used initially in Ref. Q] 
to study the redistribution phenomenon and further in 
Ref. jllj to establish its relation with the classical mon- 
odromy. We use this Hamiltonian ([[]) to illustrate the 
strict formulation of our result (j^), but its validity ex- 
tends to a general H(S,3/j). 

In the two extremes limits t = (no "spin-orbit" 
coupling), and t = 1 ("spin-orbit" coupling), the en- 
ergy level spectrum of ([!]) shows different patterns of 
energy levels into bands indexed by g. For t = 
all energies E g = g, g S {— s, ..,+s} appear with the 
same multiplicities N g = (2j + 1). For t — 1 the 
spectrum is split into degenerate multiplets character- 
ized by different eigenvalues of the coupled angular mo- 
mentum N 2 = (J + S) 2 . As in the case of standard 
spin-orbit coupling with j > s there are 2s + 1 differ- 
ent levels E g = [n{n + 1) - j(j + 1) - s(s + l)]/(2j), 
g = n — j G {— s, . . . , +s}, with different multiplicities 
N g = (2j + 1) + 2g. The two different limiting cases for 
the structures of the (2s + 1) bands suggest to introduce 
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a new quantum number C g associated with the value of 
N g within each band. 

To define C g for a general Hamiltonian H(S,3/j) we 
assume that j 3> s, so that it is physically reasonable 
to consider 3 c i = J as classical, whereas S remains 
quantum. The classical dynamics for 3 c i can be de- 
fined through the SU(2) coherent states |J ci ) The 
classical phase space for 3 c i is the sphere Sj . From 
dJ c i/dt — dj cl H c i x J c i and because of the factor 1/j 
in Eq. ([I]), J c / corresponds to a slow dynamical variable 
compared to S, and the Born-Oppcnhcimcr approxima- 
tion suggests to consider for each 3 c i, the Hermitian op- 
erator H s (3 c i) = (J c i\H\J c i) acting on H s , with spec- 
trum H s (J cl )\ip gi j cl ) = EgjJipgjJ, g € {-s,..,+s}. 
We suppose that for each J c /, the (2s + 1) eigenvalues 
are isolated: £- s ,j cl < S_ s+ i,j eI < ... < E +S< j cl . This 
is the generic situation, because degeneracies are of codi- 
mension 3, and 3 c i e Sj is only two-dimensional. 

For each level g, let us note [|^ 9 ,j c! )] the eigenvector 
defined up to a multiplication by a phase e ta . The appli- 
cation J c i — ► [IV'g.Jcj)] defines then a U(l) fiber bundle 
over the sphere Sj, which is the set of the all possible 
phases a for every values of J c ;. The topology of this 
bundle is characterized by a Chern number C g € % [[l3| . 
C g reveals the possible global twist of the fiber of phases 
a over the sphere Sj, in the same way the well known 
Mobius strip is the real line fiber bundle over the circle 
S 1 with a global twist +1. 





FIG. 1. The Chern index C g expresses the twist made by 
the circles for the phases a of the eigen-vectors e ta |^ s ,j <:; } 
over the sphere of spin J c ;, in the same way the Mobius strip 
has a twist made by the lines over a circle. 



Note that 



+ s 



(2) 



g=-s 



just because of the additivity of Chern indices, and be- 
cause the eigenvectors (l^g.Jci)) s P an the space Ti. s which 
does not depend on J c /. 

The Born Oppenheimer approach tells that the total 
spectrum H\<j>i) = E t \(f) t ),i G {1, .., (2s + l)(2j + 1)} in 
HlQ^, can be represented as formed in 2s + 1 groups of 
levels (bands) with eigen-functions of each group span- 
ning vector spaces: L g C Htoti 9 e {~ s > •••■> + s }- We 
precise now the number of levels in each group: 



Theorem. 

For a general Hamiltonian H(S,3/j), and for j large 



enough then dimi g = dim?^ — C gi i 

N g = (2j + 1) - Cg. 



(3) 



The number N g of states \<f>i) in each band L g of Ti-t i ^ s 
given by the quantum number j and an additional quan- 
tum number C g , namely the topological Chern index. 

A few remarks are in order here. Since (2j + 1) is the 
number of quanta in the classical phase space Sj for spin 
Jci, Eq.(^) looks like a Weyl formula with a correction. 
The index C g has been defined and can be computed in 
"a semi-quantal" approach where J c / is considered as a 
classical variable and S quantum. Nevertheless Eq.(||) 
provides an information on the full quantum problem: 
the spectrum of H . Finally the topological nature of C g 
reveals a qualitative and robust property of the spectrum 
of H, stable under perturbations, provided no degeneracy 
appears between consecutive bands. One can say that C g 
expresses a topological coupling between the dynamical 
variables J and S. 

Two approaches to the computation of the Chern in- 
dices of different bands will be suggested. The first one 
we use below is algebraic. The second one uses the 
Berry's connection, and is based on the curvature for- 
mula [jl3| . This last formula could be useful for numerical 
computations. 

The algebraic calculation is based on the geometric 
interpretation of the Chern index C g as the total in- 
tersection number between the one-dimensional curve 
[|^fl,Je()]jci in P{fts) with the hyperplane (IV'o))^ = 
{|<^}such that (<p\ipo) = 0}. Here |^> ) € Tt s is arbitrary. 
Each intersection has number a = +1 ( — 1) if the curve 
orientation is compatible (incompatible) with the orien- 
tation of the hyperplane. C g is the sum of these inter- 
section numbers |14]. 

The application of this algebraic method of calculation 
to the two limiting cases of Hamiltonian ([!]) is immediate. 
For t = the Hamiltonian H s (3 c i) = S z does not depend 
on 3 e i, the application J c ; — » [l^g.j^)] is constant, the 
topology of the bundle is trivial with zero Chern index: 

H). For t = 1 we have 



C g = 0, in accordance with Eq 



H S (J C 
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t(J c / ■ S) 



with 3 c i a vector on the sphere Sj, and S a vectorial op- 
erator in H.s- The eigenvector \i/j gi j ol ) is easily obtained 
from the eigenvector \m s = g) cx S'^ s+9 ''|0) of S z by a 
rotation which transforms the z axis to the 3 c i axis on 
the sphere. It is convenient to choose a coherent state 
|Soci) as t ne reference state l^o); where S c \ is an ar- 
bitrary classical spin S. An intersection of the curve 
[|V , s,J c i)]jci w ith the hyperplane (l^o)) is then given by 
the equation (So )C z|V ; fl,Je() = which has a very simple 
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interpretation: the point So. c z is a zero of the Husimi rep- 
resentation Hus(S c ;) = |(S c /|^ 9 ,j ci )| 2 of the state \i/} g ,j cl ) 
on the sphere (the classical phase space of S), shown 
in figure |[ The Husimi representation has (s — g) zeros 
inside the (oriented) trajectory, and (s + g) zeros outside 
[p~5|| - When the axis J c ; moves with direct orientation 
on the whole sphere, the (s — g) zeros pass with positive 
orientation over the fixed point So, c z, giving u\ = s — g, 
whereas the s + g zeros pass with negative orientation 
over S ,d, giving <7 2 = ~{s + g). So C g = o\ + a 2 = -2g 
in accordance with Eq. (^) . 



A J axis 




(s-g) zeros 



Classical support of state ly j > 
(= classical trajectory) 



FIG. 2. The Husimi distribution for the state |^ g ,j ai ) on 
the sphere Sf for spin S. So, c ; is a reference point. 

In this last example, the variation of Chern indices 
AC g = — 2g occurs at t = 1/2, with a degeneracy be- 
tween the bands at J c ; = (0,0, —j), giving H s (3 c i) = 0. 
In the case of two bands (s = 1/2, g = ±1/2) then 
AC g= ± = =pl, and the two bands have a conical con- 
tact at the degeneracy. In the vicinity of the degeneracy, 
and for j — ► oo, we observe that [J_, J + ] = — 2J Z ~ 2j, 
so a = J-/^/2j, a + = J + /^/2j fulfilled the harmonic os- 
cillator commutation relations. In the basis |±) = \m s = 
±1/2) of H s , the expression of H can be simplified and 
gives 



H = 



1 



—t a 
a+ t 



(4) 



(2t — l)\/2] ■ We also scales the energy with 



with t 

E = Ey/2], and note \n) oc a +n \0). The stationary 
equation H\<j>) = 
n = l,2,..., |0±) 

E± = ±v / rT 



E\(/)} can easily be solved, giving for 
_ = \ n )\-) + ^i/(E±+t)\n-l)\+) with 
t 2 , and one single state \<po) = |0)|— ) with 



Eq = t. Figure [| shows this spectrum with the simplified 
expressions of \4>n) obtained for large \t\. These simpli- 
fied expressions involving the states \n) of the harmonic 
oscillator express the quantized modes for small oscilla- 
tions near the extrema of the two bands [^6) . We clearly 
observe the exchange of one state in the spectrum, giv- 



ing AN g= ± = ±1. As a consequence, AN g + AC g 



0. 



This gives for each band a conserved quantity, namely 
(N g + Cg). The variation AC + = —1 can also be consid- 
ered as "the topological charge" associated to the degen- 
eracy of the 2x2 matrix (Q) where a c ; = x c i — ip c i is a 
classical variable. [AC+ = — 1 is also obtained from Eq. 



(^), for a little sphere in space (x c ;,p c /, t) centered at the 
degeneracy (0, 0, 0)] . 
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FIG. 3. Spectrum of rotational states for the Hamiltonian 
([|). t = corresponds a conical contact (degeneracy) between 
the two bands, with a topological charge AC+ = — 1, and an 
exchange of one rotational state, AiV+ = +1. 

For a general Hamiltonian H(S,J/j) the very simple 
model Eq. (|]) provides the general mechanism for the ex- 
change of one state in the vicinity of every degeneracy be- 
tween two consecutive bands. In the trivial case Hq = S z , 
we have computed the value of the conserved quantity: 
Ng + Cg = 2j + 1 . We deduce that it is still correct when 
H is deformed to H(S,3/j), proving formula (||). 

The second example, addresses the Chern indices of 
the two-state model corresponding to s = 1/2, i.e. to 
dimH s = 2. In a given fixed basis , say \m s — ±1/2), 
the matrix of H s (J c i) has the form: 



11 = ( fell ( Jc ') ^12 ( J c/) 

<Yl2(Jc/) h 2 2(Jcl) 



(5) 



This matrix will give two vibrational sub-bands with 
Chern indices C_ and C+. In Ref. |17| ], it is shown 
that the Chern indices have the following property: 
Let J* be the zeros of hi2(3 c i), and <r(J*), their de- 
gree defined as follows: take a small direct circle 
around J*; its image by /ii2(J c i) is a closed curve 
around with <r(J*) turns. Define the set S + = 
that (h 22 (3ci) - /Hi(Jcl)) > 0}. Then 



{J 



d € Sj such 



C 4 



E < 

j*es+ 



-C- 



(G) 



A direct consequence of this property is that a change 
of Chern index can only occur when simultaneously 
foi2(J*) = ^22(J*) - h\x{3*) = 0. This corresponds to 
a degeneracy in the spectrum of H Sl and a conical con- 
tact of the two bands. 

Formula (|^) can be applied to a Hamiltonian describ- 
ing the rotational structure of the doubly degenerated 
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vibrational state of a tetrahedral (or octahedral) spher- 
ical top molecule [jl8| . The most general Hamiltonian, 
taken up to the third degree in J, has the form (|J) with 

MJ) = (■£ - Jy) if + lXJ x J y J z /j\ 

h 22 (J)-h 11 (J) = (3J 2 z -j(j + l))/j 2 , 

with parameter X £ JR. In this case the set S + includes 
all points around north and south pole for which J z > 
TjVWTT) or J z < ^j( J + l). As h 12 (3) = for 
simultaneously J x = ±J y , and J x J y J z = there are 
two points J* G iS> + : the north pole J z — j and the 
south pole J z = —j. We consider J going through a 
closed path surrounding each J* to calculate c(J*).This 
gives for north and for the south poles a nort h = a south = 
2sign(Y), so that C+ = 4sign(Y) and C_ = -C+. 

This calculation explains why the rotational structure 
of doubly degenerate vibrational state is generally split 
into two sub-bands with respectively 2j'+5 and 2j— 3 lev- 
els [fnljsl . In our current approach the appearance of two 
bands with Chern indices ±4 for the Hamiltonian (|J) is 
due to the formation of eight degeneracies [equivalent by 
symmetry] between the two vibrational bands at X = 
parameter value. More generally the characterization of 
the rovibrational structure of molecules and its possible 
modification under the variation of some physical param- 
eters like total angular momentum can be done system- 
atically by using Chern indices as topological quantum 
numbers. In Ref. |Q similar effect was discussed without 
explicit introduction of Chern indices. 

The molecular application studied in this Letter was 
mainly inspired by the role played by topological Chern 
indices in the integer quantum Hall effect p0[ |. In this 
context, Chern indices describe the topology of Floquet 
bands (k) — » [|^ g (k))] where k is the Bloch wave vector, 
and give a quantum Hall conductance a g = (e 2 /h)C g un- 
der the hypothesis of adiabatic motion of k when a weak 
electrical field is applied. Contrary to Eq.(||), their sum 
for a given Landau level is ^ g C g = +1, because of the 
non trivial topology of the quantum space j2Q]. Many 
properties of the Chern indices are similar: a change 
AC g = ±1 is related to a conical degeneracy between 
consecutive bands. The application of semi-classical cal- 
culations of C g done for the Hall conductance in Ref. [Bl| , 
will be the subject of future work. 

In summary we have discussed the role of Chern quan- 
tum numbers to molecular spectroscopy. The interpreta- 
tion of good integer quantum numbers associated with ro- 
tational structure of different vibrational bands in terms 
of topological Chern numbers has naturally a wide appli- 
cability. These indices can be introduced any time when 
the adiabatic separation of variables enables one to split 
the global structure into bands associated with the "fast 
motion" and their internal structure described by a "slow 
motion" on a compact phase space. We have considered 
here only the problem when the dimension of the classical 



phase space formed by the slow variable 3 c i is 2. Only the 
first Chern class appears in this case. Extension to higher 
dimension requires more delicate physical interpretation 
and more sophisticated mathematical tools in relation 
with the index theorem of Atiyah-Singer. In molecu- 
lar spectroscopy many problems with intra-molecular dy- 
namics are known in great detail. They can be used to 
test the applicability of this new concept. 

We gratefully acknowledge Y. Colin de Verdiere for 
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comments. 
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